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Abstract In this article, we report the recent progress on decoherence dynamics of 
electrons in quantum dot quantum computing systems using the exact master equa- 
tion we derived recently based on the Feynman- Vernon influence functional approach. 
The exact master equation is valid for general nanostructure systems coupled to multi- 
reservoirs with arbitrary spectral densities, temperatures and biases. We take the dou- 
ble quantum dot charge qubit system as a specific example, and discuss in details the 
decoherence dynamics of the charge qubit under coherence controls. The decoherence 
dynamics risen from the entanglement between the system and the environment is 
mainly non-Markovian. We further discuss the decoherence of the double-dot charge 
qubit induced by quantum point contact (QPC) measurement where the master equa- 
tion is re-derived using the Keldysh non-equilibrium Green function technique due to 
the non-linear coupling between the charge qubit and the QPC. The non-Markovian 
decoherence dynamics in the measurement processes is extensively discussed as well. 
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1 Introduction 

The investigation of quantum coherence and quantum entanglement in quantum dot 
systems has been attracting much attention in the past decade because of its scalabihty 
for quantum computer. Due to the rapid development of nanotechnology, implement- 
ing quantum information processing using nanostructures is very prospective. Among a 
diversity of nanostructures, a prototypical one is made by fabrications on the interface 
of a heterostructure semiconductor as a gate-defined region containing a set of coupled 
electronic states. Electrodes are implanted around this region to control the bias across 
it, and also several gates are implemented to adjust the electronic states within the 
central area as well as their couplings to the surrounding electrodes. Because the elec- 
tronic structures of these systems can be well controlled, they have been intensively and 
extensively studied for many purposes including understanding fundamental physical 
issues and building useful quantum devices. However, it is also because of its openness 
to the controlling electrodes, the external gates and various degrees of freedom in the 
host material, the central electronic region is intimately entangled with the surrounding 
environment. Quantum entanglement among the basic units for quantum information 
processing, i.e. quantum bits or simply qubits, is the unique resource for exponentially 
speedup over classical algorithms. However, quantum entanglement between the qubit 
system and its environment is indeed a debt to quantum information processing. The 
entanglement between the system and its environment induces severe loss of quantum 
coherence of qubits and destroys eventually the useful entanglement among qubits. 
It has been realized that decoherence is the most difficult obstacle for any successful 
quantum information processing. Therefore, the investigation of decoherence dynamics 
due to the entanglement of the qubit systems with the surrounding environment and 
the development of decoherence control protocols become the central issue in the study 
of quantum information processing. In this article, we will present a general theory for 
the description of the real time dynamics of a central electronic system coupled to 
multi-electrodes. The theory is build on the exact master equation for general nanos- 
tructures we derived recently using the Feynman- Vernon influence functional approach 
[l]. This exact master equation is capable to depict the detailed decoherence dynamics 
of electronic states in the nanostructure with arbitrary spectral density at an arbitrary 
initial temperature of the reservoirs and an arbitrary bias applied to the reservoirs. 

Historically, since it was first proposed by Feynman and Vernon in 1963 [5] for 
quantum Brownian motion (QBM) modelled as a central harmonic oscillator linearly 
coupled to a set of harmonic oscillators simulating the thermal bath, the influence 
functional approach has been widely used to study dissipation dynamics in quantum 
tunneling problems [3] and decoherence problems in quantum measurement theory 
[31[S]. In the early applications, the master equation was derived for some particular 
class of Ohmic (white-noise) environment [6]. The exact master equation for the QBM 
with a general spectral density (color-noise environments) at arbitrary temperature 
was obtained by Hu and co-workers in 1992 [7]. Applications of the QBM exact mas- 
ter equation cover various topics, such as quantum decoherence, quantum-to-classical 
transition, quantum measurement theory, and quantum gravity and quantum cosmol- 
ogy, etc. [8ll9lll0j. Very recently, such an exact master equation is further extended to 
the systems of two entangled harmonic oscillators [TT] and two entangled optical flelds 
[12lll3j for the study of non-Markovian entanglement dynamics. Nevertheless, using 
the influence functional approach to obtain the exact master equation has been largely 
focused on the bosonic environments in the past half century. The extension of the 
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influence functional approach to fermion environments is just started [T], where the 
Feynman path integral in terms of fermion coherent states [14] must be used. 

In fact, the investigation of the environmental efl'ects on the subsequent dynamics 
of the principal system one concerned is one of the essential issues in many currently 
interesting research topics related to quantum decoherence and dissipation phenomena. 
The quantum decoherence and dissipation phenomena are also associated with quan- 
tum non-equilibrium dynamics, a subject that has been investigated over a half century 
yet is still not completely understood. Quantum decoherence in quantum computa- 
tions and quantum information processes is one of the biggest challenges for practical 
applications of quantum mechanics. Using the Feynman- Vernon influence functional 
approach to understand decoherence phenomena is fundamental in the sense that it 
allows to completely integrate out the environmental degrees of freedom. As a result, 
the back-action of the environment to the system can be fully taken into account. The 
resulting non-perturbation (with respect to the coupling between the system and its 
environment) master equation is indeed desirable for quantum information processing 
since the fast manipulations to qubit states requires a strong coupling among the con- 
stitutes. It enables us to precisely explore the real-time dynamics of the electron charge 
coherence in various different nanostructures under different manipulating conditions. 

On the other hand, quantum measurement (readout of qubit states) also often 
induces severe decoherence. Therefore, precision readout of qubit states is another 
big challenge in quantum information processing. Again due to the rapid progress in 
nanotechnology, electron state readout in nanostructures has been investigated exper- 
imentally with qubit systems coupling to a mesoscopic measurement device, such as 
single electron transistor (SET) or quantum point contact (QPC). Not only its prac- 
tical application to quantum information processing, but also the theoretical interests 
in the measurement-induced quantum decoherence have attracted much attention. In 
particular, in the investigation of quantum dot quantum computing with charge qubits, 
the QPC has been served as an ultrasensitive electrometer. The theoretical understand- 
ing of charge qubit measurement through the QPC was previously treated based on 
the Markov approximation [15lll6lfT7] in which the time scale of the qubit dynamics 
is assumed much longer than that of the tunneling-electron correlation in the QPC. 
The treatments based on the Markov approximation, in fact, only describe the qubit 
dynamics in the time-asymptotic quasi-equilibrium state. However, for the quantum 
measurement in terms of SET or QPC, the qubit decoherence occurs in the time scale 
of the same order of the tunneling-electron correlation time in the measurement equip- 
ments, where the non-Mar kovian dynamics of the qubit is significant. In the second 
part of this article we shall discuss a fully non-equilibrium theory we have developed 
recently |18II19| for understanding the qubit decoherence induced by the QPC mea- 
surement, based on the Keldysh non-equilibrium Green function technique |20) . 

The Keldysh non-equilibrium Green function approach is well developed for a per- 
turbation treatment to the non-equilibrium dynamics in many-body systems [21ll22j . 
This approach has been used to study various transport phenomena in nanostructures 
[23j . To make the non-equilibrium effect of the electrical reservoir more transparent, 
a real-time diagrammatic technique was constructed to diagrammatically calculate 
correlation functions of the electrical reservoir order by order in the perturbation ex- 
pansion [24ll25j . In the investigation of the non-Markovian dynamics in the double dot 
charge qubit induced by the QPC measurement, we developed an alternative real-time 
diagrammatic technique [18]. The master equation for the charge qubit dynamics is 
derived and expressed in terms of irreducible diagrams up to all orders. The effect of 
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the fluctuant measurement reservoir due to the interaction with the qubit system can 
be fully taken into account. The measurement-induced non-Markovian decoherence in 
the qubit dynamics can be explicitly studied in this formulism. 



2 Exact master equation of electron transport systems 

The prototypical nanostructure under consideration can be generally described by the 
following electronic Hamiltonian 

ij ij 

+ J2^kak4cak + t*akclkdi) , (1) 
ak iak 

where the first two terms are the general Hamiltonian Hd of electrons in the central 
region of the nanostructure. For simplicity, we shall ignore the electron-electron inter- 
action in this article, i.e., let Uij = 0. The third term is the Hamiltonian Hj^ describing 
the noninteracting electron reservoirs (the source and drain electrodes, etc.) labeled by 
the index a. The last term is the electron tunneling Hamiltonian Hj^ between the 
reservoirs and the central system, dj(di) and c^^^j^icak) a-re the electron creation (anni- 
hilation) operators of the central system and the surrounding reservoirs, respectively. 
Also, throughout this work, we set h — 1. The initial inverse temperature for the a-lead 
is denoted as Pa = l/k^Ta- 

Using the Feynman- Vernon influence functional approach extended to fermion co- 
herent state representation [14] . we can integrate out completely the degrees of freedom 
of the electron reservoirs and obtain the following exact master equation [1], 

p{t) = -^[HDit),p{t)] 

ij 

+ r['^{t)[d,p{t)di ~ d\p(t)d, ~ d\d,p{t) + pmAVi ■ (2) 

This master equation determines completely the quantum coherence dynamics of the 
central electronic system. p(t) is the reduced density matrix of the central system at 
time t after the environmental degrees of freedom are completely eliminated. H'j-){t) — 
"Ylij ^ij{^)4dj is corresponding renormalized Hamiltonian. Other two non-unitary 
terms in the master equation describe the reservoir-induced dissipative and noise pro- 
cesses with time-dependent dissipation and fluctuation coefficients, r{t) and r^{t). All 
these time-dependent coefficients in the master equation are derived non-perturbatively 
and are exact. The time-dependent fluctuations of the energy levels, ^'ait), and the tran- 
sition couplings between states, e^j(t) for i ^ j, are the renormalization effects that 
contain the energy level shifts and the coupling changes between them, while r{t) and 
r^{t) are the dissipation and fluctuation effects that depict the full non-Markovian 
decoherence dynamics. These are all the back-action effects risen from the electron 
tunneling processes between the central system and the reservoirs. 
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The explicit time dependencies of these coefficients in the master equation are given 
as follows 

e'^jit) = ^[u{t)u-\t) - {u^r\t)u\t)Uj, (3a) 

nj{t) = -^[u{t)u-\t) + {u^r\t)u\t)],j, (3b) 

r|(t) = [u{t)u-\t)v{t) + v{t){uY\t)u\t) - «(t)],,. (3c) 

The elementary functions u{t) and v{t) are time-dependent matrices following the 
dissipation-fluctuation integral-differential equations 

u{T) + ieu{T) + J2 [ dr'Fa{r -r')u{T') = 0, (4a) 

v(T)+iev{T) + J2 f drFair -t')v{t') = J2 f dr (r - t')u{t'), (4b) 

subjected to the boundary conditions w(fo) = I a-nd w(fo) ~ 0, where ■u(r) = u^{t — 
r -\- to). Introduce the spectral density Jj"(e) ~ "^^i^^attak^jakSi^ — ^ak)) then 
the integration kernels in the above equations can be expressed as Fijij — r') = 

/ #J5 (.)e— and i^^(r-r') = / |^/a(e)J^-(.)e— where /„(.) = 
— o— ; — - — r— r is the initial fermi distribution function of the reservoir a held at the chem- 
ical potential /Xa- These are the so-called dissipation-fluctuation kernels and represent 
dissipation and fluctuation effects on the central system due to its coupling to electron 
reservoirs. The relation between these two kernels are then seen from the dissipation- 
fluctuation theorem. The real time non-equilibrium electron dynamics of the central 
systems is fully described by the master equation. 

3 Coherence control of double dot charge qubit 

In the quantum computing scheme in terms of double quantum dots where the electron 
charge degree of freedom is exploited, the effects in deviating the coherency of charge 
dynamics are usually summarized in the fluctuations of the inter-dot coupling and en- 
ergy splitting between the two local charge states. The amplitudes of these fluctuations 
may be estimated from measurements of the noise spectrum of electron currents and 
the minimum line width of elastic current peak [26]. Parallel theoretical works have 
been developed with different approaches in the literature for the purposes of both 
simulating the experimental results and understanding the physical mechanisms living 
in the double-dot. With the master equation shown above, a complete investigation 
can be carried out for the decoherence dynamics of the double dot charge qubit under 
coherence controls. The decoherence of the electron charges can be addressed with the 
back-reaction effects of the reservoirs to the double dot (including not only the fluc- 
tuations of the inter-dot coupling and energy splitting between the two local charge 
states but also the non-unitary dynamical effect induced by dissipations and noises) 
being fully taken into account. 

For the double dot charge qubit (a schematic plot see Fig. [TJ , the corresponding 
Hamiltonian of the central system becomes 

Hdot = Eid\di + E2d\d2 -f nQ{d\di + d\d2) , (5) 
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Charge Qubit 




Fig. 1 (color online). A schematic plot of the double quantum dot coupled to the source and 
drain electrodes. 



where Ei and E2 are the on-site energies of the left and the right dot respectively, 
J7o is the tunnel coupling between the two sites. The double quantum dot is coupled 
to two electron leads on the left and the right side. Specifically, the left dot is only 
coupled to the left lead while the right dot is only coupled to the right lead. To closely 
monitor the time evolution of the electronic states in this system, it is more convenient 
to specify the master equation in the particle number basis of the no electron state 
|0), the state of one electron in the left dot one electron in the right dot |2), and 
the double occupied state of two electrons in two dots |3). Then the master equation 
is reduced to a set of coupled rate equations in this basis. Furthermore, the double dot 
charge qubit works in the single-electron regime, where the strong Coulomb blockade 
(Coulomb repulsion between the two sites) plays an important role. In other words, 
the doubly occupied state |3) must be excluded physically. This can be easily taken 
into account by excluding the double occupancy from the rate equations properly. The 
consequent rate equations excluding the double occupancy read [T] 

POO = AlPll + r2ipi2 + A2/521 + r22P22 + Trf^poo, (6a) 
Pll = -Alpll + -S"* Pl2 + S-P21 - rf^Poo, (6b) 

P22 = -r22p22 + S'+P12 + S+p21 - r^2P00i (6c) 

P12 = [-je' - ^trrlpi2 + S'+pii + S-P22 - r^2P00, (6d) 

where f(t) = 2r{t) + F^^ (t), H±(t) = ±iei2(t) - ^A2(i) and e'{t) = e'n(t) - €'22(1), 
r{t) and r^{t) are the time-dependent coefHcient matrices in the master equation ((2| 
with i = 1,2. We can reproduce other rate equations from (O in the corresponding 
Born-Markov limit used in the literature [27l[28][29] . 

To study the decoherence dynamics of the charge qubit, the spectral density must 
be specified. Here we use a Lorentzian spectral density. For the double dot charge qubit 
system, since the left (right) dot only couples to the left (right) electrode, the spectral 
density can be written explicitly as 

^^_j,^^2 + wi^V, (7) 

with a — L{R) for i — 1(2). The Lorentzian spectral widths W^^fj axe indeed the 
bandwidths of the densities of states for the source and drain, respectively. They depict 
the characteristic times of the corresponding reservoirs. The larger the widths Wi^ n 
are, the shorter the correlation time of the reservoirs give and the lesser the memory 
effect shows. In the wide band limit (WBL) that is widely used in the literature, the 
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Fig. 2 (color online). The time dependence of transport coefficients in the master equation 
where A'{t) = 2e'j^2(*)- The inter-dot tunnel coupling Qq is set to be 5/ieV with the fixed 
chemical potentials at fj, — E = 50/ieV and the Lorentzian spectral width W = 20fieV . The 
red line, the blue dashed line and the black dot lines correspond to 7^ = 5, iO and 25^eV, 
respectively. 



spectral density becomes energy independent. As a matter of fact, all the memory effect 
is wiped out. The couphng constants Fi^ jj are the tunneUng rates of electron tunneled 
into the double dot from the source and the drain, respectively. The tunneling rates 
describe the strength of the tunneling process between the dot and the reservoirs 
covering the leakage effect, and are experimentally adjustable through gate voltages. 
Varying the coupling strength to the environment and the correlation time scale of the 
reservoirs, we are able to analyze realistically the decoherence dynamics of the double 
dot charge qubit. 

For manipulating the charge qubit coherence, the device is set at the resonant con- 
dition El = E2 = E. The fermi surfaces of the two electrodes are aligned = /x/j = /x 
above the resonant level ^ — _E > so that electrons are kept in the dot with minor 
probability of leakage. To simplify the problem, we consider only the symmetric dou- 
ble dot system where F]^ = = F and = Wfi = W. The electron temperature 
is kept at 100mA' throughout the rest of the calculations. The time-dependence of 
the transport coefficients in the master equation is essential for the non-Martkovian 
dynamics during the coherence controls of charge qubit states. Here we first examine 
the time dependent behavior of these transport coefficients. Fig. [2] shows an example 
of the time dependencies of those coefficients at a few different strengths of coupling 
(tunneling rate) between the qubit and electron reservoirs, where the coupling strength 
can be controlled by tuning the gate voltages in experiments. The time is in unit of 

To — 27r/y^ {Ei — i?2)^ + 4:^q = tv/Hq. As one can see the stronger the coupling be- 
tween the qubit and the environment is, the larger fluctuation the coefficients shows 
that reveals the non-Markov feature. As it is known the rate between the bandwidth W 
and the tunneling rate F tells how important the non-Markovian dynamics can be. The 
width W characterizes the correlation time of the electron reservoirs, a large width W 
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Fig. 3 (color online). The time evolution of the reduced density matrix elements for the charge 
qubit solved from the master equation. All the parameters are the same as that given in Fig. |2] 



corresponds to a short characteristic time of the reservoir. Therefore, a smaller W/F 
will make the non-Markovian process more obvious in electron dynamics. This fea- 
ture is shown by the time-dependence of the transport coefficients, as we can see from 
Fig. [21 When W/ F becomes large enough, the time-dependence of the transport coeffi- 
cients will disappear and the Markovian limit may be reached. Equivalently speaking, 
a larger width W corresponds to a shorter correlation time of the electron reservoir, 
the coefficients at larger width reach their steady limit sooner. 

Correspondingly, the dynamics of the charge qubit is also examined at different cou- 
pling strengths by plotting the time evolution of the reduced density matrix elements. 
Fig.[3ja)-(d) show the time evolution of the electron population occupying the first dot 
given by pii(t), the probability of empty dots, i.e. the leakage effect by poo(i)) ^nd the 
real and imaginary parts of the off-diagonal matrix element pi2{t). For the symmetric 
dots {El = E2) concerned here, the initial energy eigenstates of the double dots are 
given by the bonding and anti-bonding states, |±) = i |2))- The relaxation 

and decoherence of the qubit states are then determined by the decays of the matrix 
elements = ^(pii -I-P22) +Repi2 and {-l-|p|-> = ^(pii - P22) - jlmpi2, respec- 

tively. In other words, the decay rate of Repi2(t) is the relaxation time Ti (because 
Pll + P22 — 1) while the decay rates of Pii(t),P22(i) or Impi2(t) give the decoher- 
ence time F2, as we have discussed in [1]. These properties are clearly shown in Fig.[3l 
where the oscillation decays in Fig.[3ja) and (d) depict the decoherence dynamics while 
Fig. EJc) describes the relaxation process (energy dissipation) and Fig. [Sfb) shows the 
small leakage effect. We find that the decay of the charge coherent oscillation is well 
described by a simple exponential decay for the off-diagonal reduced density matrix 
elements. The diagonal matrix elements (populations) are better described by a sub- 
exponential law when charge leakage is not negligible, otherwise simple exponential 
decay is better for all the cases. The relaxation time Ti and the decoherence time T2 
can be extracted from the exact numerical solution Repi2(t) and Impi2(t), respectively. 



9 



Tl „=20fieV, Wl „=25neV „=1 OneV, „=25neV 




01234 01234 

t/T„ t/r„ 




Fig. 4 (color online). A comparison of the exact solution (red lines) with the Markov limit 
(blue dashed linos), where the inter-dot tunnel coupling J7o = 5/ioV with the fixed chemical 
potential at fj, ~ E = 50^cV. 



with the result Ti < 2T2 being of the order of a few nanoseconds or less for a broad 
parameter range we used pp. 

Furthermore, a stronger tunneling coupling F to the electron reservoirs (compare 
to the Lorentzian spectral width W) leads to a faster decay of the coherent charge 
qubit oscillation and also causes severer shifts in the charge oscillation frequency. The 
non-Markovian decoherence dynamics of charge qubit is dominated by two major ef- 
fects, the memory effect and the leakage effect in the double-dot gated by electrode 
reservoirs. The former becomes a dominate effect when the time scale (the inverse of 
the width W) of the reservoirs is comparable to the time scale (~ Tq) of the double- 
dot. The latter becomes important when the electron tunneling strength between the 
reservoirs and dots is tuned to be large. Strengthening the couplings between the reser- 
voirs and dots disturbs the charge coherence in the double-dot significantly. However, 
reasonably raising up the chemical potentials fiL,fl. can suppress charge leakage and 
maintain charge coherence. The left uncontrollable decoherence factor is the spectral 
width which characters how many electron states in the reservoirs effectively involving 
in the tunneling processes between the dots and the reservoirs. It is quite interest to 
see that in any case, the Markov limit (based on the rate equations derived by Bran- 
des and Vorrath under Born-Markov approximation but without including the phonon 
coupling [29]) deviates largely from the exact solutions, as shown in Fig.|4l It indicates 
that quantum information processing is indeed non-Markovian with the double dot 
charge qubit. 



4 Charge qubit with QPC measurement 

In this section, we present a study of the charge qubit measurement using QPC in 
the tunnel junction regime. The transmissions of all tunneling channels cross the QPC 
barrier are required to be small enough in this regime such that the electron tunneling 
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Fig. 5 A schematic plot for the QPC measurement of double dot charge qubit. 



becomes sensitive to the qubit state. The qubit information can then be extracted from 
the output signal of the QPC, while the back-action of the measurement to the qubit 
states can also be taken into account. In order to explore the qubit decoherence induced 
by the QPC measurement, the corresponding Hamiltonian takes the form [15] 

H = H^ot + ^eaka^^kaak + ^{<lkk' a^Rk-'^Lk + qlk'°''^Lk°'Rk')- (8) 
a,k kk' 

where -ffdot is the Hamiltonian of the double dot charge qubit given by Eq. The 
second term is the Hamiltonian of the QPC electronic reservoirs consisting of the source 
(a = L) and the drain (a = R) electrodes with the energy levels e^.^,, and the chem- 
ical potentials ii^ a^nd /i/j, respectively. The last term is the interaction Hamiltonian 
describing the electron tunneling processes through the QPC with a qubit-state depen- 
dent hopping amplitude qkk' = ^kk' ~ SQkk'd\di. The single electron occupies on the 
first dot leads to a variation in the barrier of the QPC, and the electron hoping am- 
plitude is thus modified from flkk' to 0'f.f., = flkk' ~ Sfikk' ■ The charge qubit state is 
measured through the electron tunneling across the source and the drain. A schematic 
device is plotted in Fig. [5l 

To explore the measurement-induced decoherence of charge qubit, the QPC can be 
considered as an environment coupled to the charge qubit. Then the qubit dynamics is 
determined by the master equation for the reduced density operator p{t) =trjf [ptot{t)], 
where ptot (t) is the total density operator for the whole system of the qubit plus QPC 
reservoir, and the partial trace tr^ integrates over all the degrees of freedom of the 
QPC reservoir. Since the qubit-reservoir coupling involves two-body interactions, the 
influence functional approach for integrating over all the degrees of freedom of the 
QPC reservoir cannot be exactly carried out. Therefore, we developed a real-time dia- 
grammatic technique based on the Keldysh non-equilibrium Green function approach 
to study the decoherence effect of non-equilibrium electrodes on the qubit. The master 
equation for the charge qubit coupled with the QPC measurement has been derived in 
terms of the irreducible diagrams to all the orders in the perturbation expansion [181 

US], 

p^t) = -^[Hs,p{t)] -J2 drK\^''\t - r) * p(r), (9) 

where the kernel expansion K^^^\t — r) * p(r) consists of all the irreducible diagrams 
containing 2n -f 2-particle correlations. The irreducible diagrams are defined by all 
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Fig. 6 Irreducible diagrams: (a) and (b) are the leading K^^ [t — r) and the second leading- 
order K^^. ' (t — r) contributions, respectively. 

connected topology-independent diagrams along the real time axis. The connected di- 
agram means that each loop in the diagram should intersect with other loops at least 
once. Physically, in the tunneling processes classified by k'>^"'^ [t — r), there must exist 
at least two particles correlating together in arbitrary time interval of the time pe- 
riod from to to t. The leading and the second-order irreducible diagrams are shown in 
Fig. [51 where each vertex depends on the qubit operator via qi^j^i and the connecting 
free propagator specifies the non-equilibrium effect of the tunneling electron QPC on 
the qubit through the Keldysh Green function matrix. For the detailed diagrammatic 
rules and the classification of topology-independent irreducible diagrams, please refer 
to [TSlfTg]. 

Since the QPC is often used for continuous measurement of charge qubit states, 
the coupling between the qubit and the QPC can be treated as a weak coupling. For 
the study of this QPC measurement- induced non-Markovian decoherence to the qubit 
states, it may be sufficient to consider only the leading order contribution to the master 
equation. The master equation (up to the leading order) is given as follows 

p{t) = ~i[Hs,p{t)\~ f' dr[Ro,[[K{t-r),p{r)]]], (10) 
Jto 

where p{r) = e~'^-'^^*~''^ p{r)e'-"''^*-''\ the douhlehracket [[A,B]] = AB-{AB) + , the 
operator _Ro = -R(O) and R{t) is defined as 7?(t) = cos6'(|e) (e|-|5i) (pi)— sin 6'(e*'''* |gi) (e|-|- 
g-«7t (^g^^ with 7 = 4J?Q + (El — £2)^ the energy difference between the ground 

state \g) and the excited state |e) of the qubit and 6 = cos~^[(-Ei — £2)/^]- The op- 
erator K{t) = k{t)R(t), and k(t) the reservoir correlation function that characterizes 
the QPC structure. It can be expressed explicitly as 

0^e^'-''^'{fLm-fR{e')]J{e,e')+fR{e)[l-fL{e')\J{e',e)}, (11) 

where fL,R{^) = 1/(1 -I- exp /?(e — fJ.L,R)) are the Fermi-Dirac distribution functions for 
the source and the drain, J(e, e') is the spectral density of the QPC reservoir which is 
defined by 

J{e,e')=7T^gL{e)gR{e')\5n(e,e')\^ , (12) 

with gL,R{^) being the density of states of the source and the drain, and 5i7(e, e') the 
difference of the energy-dependent hopping amplitude (between the source and the 
drain) without and with the occupation of the first dot. 



12 



0.3 



25 












3 



0.0 








/ (in units : l/i2^,) 



10 



Fig. 7 (color online). The amplitude of the reservoir time correlation function during the 
measurement with different width W. Here we take the parameters w = 2Qo, P = l/^^Oi 
Vd = 20no and Fd = J7q/4, respectively. The amplitude is plotted in units SIq. 



A close connection between the qubit decoherence and the tunneling-electron fluc- 
tuation is revealed through the reservoir time correlation function k{t — r) of Eq. (Illf) . 
The reservoir time correlation function describes the variation of the tunneling-electron 
correlation during the measurement. As we mentioned in the introduction, the previ- 
ous investigations for the QPC measurement-induced decoherence of charge qubit are 
mainly focused on the Markov limit |15II16|[T7] in which the time scale of the qubit dy- 
namics is assumed much longer than that of the tunneling-electron correlation in the 
QPC. Such theoretical treatments based on the Markov approximation describe only 
the qubit dynamics in the time-asymptotic quasi-stationary state. Mesoscopically, the 
qubit decoherence occurs in the time scale of the same order of the tunneling-electron 
correlation in the QPC, where the non-Markovian dynamics of the qubit is signiflcant 
[19j . To study the qubit decoherence in the non-Markovian regime, a Lorentzian-type 
spectral density for the QPC reservoir is considered here 



where the constant specifies the decay rate of the qubit state due to the interaction of 
the qubit electron with the QPC, and the factor in Eq. (|12p has been absorbed 

into r^. The spectral width W determines the correlation time scales of tunneling 
electrons across the QPC barrier, and w is an additional parameter to characterize 
the variation of the QPC internal structure due to the interaction with the qubit 
electron. In the literature, the wide band limit, VK ^ 7 [so that J(e,e') — > Fj], has 
been often used to specify the QPC reservoir structure which simply leads to the 
Markov limit at the time-asymptotic quasi-stationary state |15II16|[T7] . In Fig. [T] the 
reservoir time correlation function is plotted. The corresponding correlation time scale 
is determined by the half width of the profile. The wider the half width of \k{t)\ 
is, the longer correlation time scale the tunneling-electron fluctuation maintains. The 
correlation time scale increases with W decreasing. 

To make the qubit decoherence feature apparent due to the QPC measurement, we 
concentrate again on the charge qubit with symmetric coupled dots Ei — £2. Thus 
the energy scale of qubit dynamics is simply characterized by 7 = 2S7o {9 = 7r/2). The 
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Fig. 8 (color online). The charge qubit dynamics for different width W . The charge qubit is 
set initially in the |1) state. The other parameters are taken as the same as in Fig. [7] (a) and 
(b) for the charge qubit dephasing, and (c) the charge qubit relaxation. 



master motion for the reduced density matrix becomes, 



peg{t) = -i'ypeg[t)-2 I dr | fc(t - r) | cos[<j!-(f - r)] {peg (r) - pge (r) } , (14a) 

J to 

Pgg{t)=2 f dT\k{t~r)\{co8[<f>~{t^T)]pee{r)-cos[<t>+(t-T)]pgg{r)}, (14b) 

J to 



where the matrix elements are defines as pijit) — {i\p{t)\j) with being the ground 
state or the excited state of the qubit, (fi{t) is a time-dependent phase of the reservoir 
correlation function k{t — t) — \k{t — T)\e~^'^^^~'^\ and we have also defined 4>^ = 
(j>{t — t) ± {t — r)7. The qubit dynamics with different W are simulated in Fig. |8] As 
we can see, with a very large width W{— 10J?o) the charge qubit undergoes severe 
decoherence, where the time correlation of tunneling-electron fiuctuation is smeared, 
and the reservoir memory effect on the qubit dynamics is washed out. In other words, 
the qubit decoherence dynamics is Markov in this regime, which almost coincides with 
the result of the Born-Markov approximation obtained in the Uterature [15lll6lfT7] . 
see the black curves in Fig. |8] However, if the QPC structure can be recasted such 
that the width W is comparable to the energy scale of the charge qubit, then non- 
Markovian processes becomes dominate in the qubit dynamics. Tunneling electrons in 
this case propagate with a larger correlation time. The charge qubit simply undergoes 
a oscillation, see Fig. |8] It is very interesting to see that the qubit decoherence is 
indeed suppressed in this case. On the other hand, when a large bias voltage is applied, 
it usually shortens the correlation time of the tunneling-electron fluctuation. A large 
amount of electrons tunneling across the QPC barrier at a large bias leads to a sharp 
profile for the time correlation function \k{t) \ such that the qubit dynamics evolves into 
a severe decoherence. However, it is interesting to see from Fig. |8] that even though a 
large bias voltage (V^i = 20i7o) is applied, the qubit decoherence can still be suppressed 
if the QPC is fabricated to reach a narrower profile of the spectral density. The qubit 
states can then be read out by a significant current signal without much decoherence. In 
conclusion, more useful non-Markovian dynamics can be manifested when the explicit 
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effect of a realistic spectral density with a finite correlation time scale is taken into 
account, where the qubit decoherence can be suppressed during the QPC measurement. 

5 Prospection 

In this article we have made an extensive discussion on the non-Markovian decoherence 
dynamics of a double dot charge qubit induced from the coherence controls and the 
QPC measurement, using the master equation we obtained from the Feynman- Vernon 
influence function approach and the Keldysh non-equilibrium Green function tech- 
nique. For quantum information processing, decoherence dynamics of individual qubit 
is no doubt a key concern in these electronic systems. On the other hand, decoherence 
dynamics of entanglement between two qubits [30ll31j is another important issue for 
quantum information processing. Some progress on decoherence entanglement dynam- 
ics in cavity optical systems has been reported based on exact master equations [11II12I 
I13j . Further applications to the decoherence dynamics of electron charge and spin en- 
tanglement in terms of the quantum dots |32II33) are expected. The full non-equilibrium 
electron dynamics and the real time monitoring of spin polarization processes in quan- 
tum dots are future topics. These together with other physical properties in various 
nanostructures, such as Kondo effect and Fano resonance, etc., as well as the transient 
dynamics of electronic quantum transports are under progress. In short, the theory we 
discussed in this work can be used to study not only the problem of decoherence but 
also many other interesting physical phenomena in various nanostructures. 
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